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We review the Random First Order Transition Theory of the glass transition, emphasizing the 
experimental tests of the theory. Many distinct phenomena are quantitatively predicted or explained 
by the theory, both above and below the glass transition temperature T g . These include: the 
viscosity catastrophe and heat capacity jump at T g , and their connection; the non-exponentiality 
of relaxations and their correlation with the fragility; dynamic heterogeneity in supercooled liquids 
owing to the mosaic structure; deviations from the Vogel-Fulcher law, connected with strings or 
fractral cooperative rearrangements; deviations from the Stokes- Einstein relation close to T g ; aging, 
and its correlation with fragility; the excess density of states at cryogenic temperatures due to two 
level tunneling systems and the Boson Peak. 
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I. INTRODUCTION 

School children have their earliest exposure to the sub- 
ject of physical chemistry when they hear about "the 
states of matter." They are taught there are gases, liq- 
uids, and solids. Van der Waals revealed to scientists that 
gases and liquids differed actually only quantitatively [1] . 
The rigidity of solids, which defines them macroscopi- 
cally, on the other hand, has been usually traced to their 
qualitatively different, periodic ordered structure [2] - an 
idea that already occurred to Kepler [3]. This idea is 
correct for crystals. The existence of structural glasses, 
i.e. amorphous substances that are rigid, calls this under- 
standing into question, however. The utility and adapt- 
ability of glasses arises from the way their properties de- 
pend on their preparation history and their seeming con- 



tinuity with the supercooled liquid state. In contrast to 
dilute gases and crystalline solids, where the properties 
can be directly inferred from the intermolecular forces, 
liquids have generally enjoyed a reputation of mystery 
among physical chemists. Until recently, the transition to 
the "glassy state" was deemed specifically one of the most 
obscure enigmas by many in the theoretical physics and 
chemistry communities [4-6]. Despite this reputation, a 
constructive molecular theory of structural glasses and 
the dynamics of supercooled liquids has been developed. 
This theory starts with the intermolecular forces. While 
it would suffer from all the well-known issues of micro- 
scopic modeling, if it were used to predict glass transition 
temperatures from scratch, for example, the theory does 
explain extremely well the known (and rather unusual) 
phenomenology of supercooled liquids and the strange 
properties of glasses at low temperatures. The theory 
also does make a large number of quantitative predic- 
tions without adjustable parameters that are borne out 
quite well by experiment. The purpose of this article is 
to briefly explain the basic ideas of this theory and how 
its predictions compare with experiment. 

The basic idea of the theory of structural glasses is to 
consider them to be aperiodic crystals. Since Bernal's 
time at least, it has been appreciated that an aperi- 
odic structure can be mechanically stable [7]. Strictly 
speaking, Kepler's assumption of periodicity would not 
be needed for rigidity, were it not for thermal motions. In 
contrast to periodic crystals, there are extremely many 
aperiodic structures and it is hard to see why they should 
have any difficulty interconverting thereby allowing flow. 
But we know they indeed do have such a difficulty since 
the flow of glass is nearly imperceptible. It is this fact 
that a theory of the glass transition must explain. We 
shall see the theory of structural glasses then generalizes 
the theory of an ordinary first order freezing transition, 
which gives a periodic crystalline structure, to a very 
much analogous change that involves freezing into a set 
of aperiodic structures which are statistically distributed 
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FIG. 1: Regimes of the aperiodic condensed molecular phase 
are shown, ranging between a dilute gas and a frozen glass. 
T v is the vaporization temperature, T m the melting point. 
Ta represents the temperature signalling the crossover to ac- 
tivated motions, which is usually but not always below T m . T g 
is the glass transition temperature which depends on the time 
scale of measurement. Below T g the system is out of equilib- 
rium and ages. Tk is the Kauzmann temperature (see text). 
Td is the Debye temperature which signals the quantization 
of vibrational motions. Below To /30, or so, the thermal prop- 
erties of the system can be phenomenologically described as 
arising from a collection of two level systems. Just above this 
point, additional quantum excitations, sometimes called the 
Boson peak, are present. 
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FIG. 2: The viscosities of several supercooled liquids are plot- 
ted as functions of the inverse temperature. Substances with 
almost- Arrhenius-like dependences are said to be strong liq- 
uids, while the visibly convex curves are described as "fragile" 
substances. The full dynamic range from about a picosecond, 
on the lower viscosity side, to 10 4 seconds or so when the vis- 
cosity reaches to 10 13 poise. This figure is taken from Ref. [6] . 



in energy. We call this change a "random first order tran- 
sition." 

We will begin this review by discussing a small number 
of key experimental signatures of the glass transition in 
Section II. In Section III, we construct the microscopic 
picture of the glassy state and the transition to it from a 
supercooled liquid, following the random first order tran- 
sition theory. A variety of temperatures characterizes 
glasses and liquids in this theory. They are graphically 
summarized in Fig.l. We will define these scales more 
precisely in the discussion below and we recommend the 
reader to often refer to this figure. Starting with a one- 
component gas, one may cool it down and compress it 
until it condenses below the critical point, T V) usually 
above the crystallization temperature T m . In this tem- 
perature range, an effective description in terms of col- 
lisional transport is valid: a liquid is just a very dense 
gas held together by an average attractive force. No two 
molecules are likely to reside near each other for any sig- 
nificant time. The time scales for molecular permuta- 
tions and collisions are comparable in this regime. All 
the pertinent information about particle-particle interac- 
tions may be encoded in low order correlation functions 
that may be computed or extracted experimentally from 
scattering experiments. In a supercooled liquid, on the 
other hand, molecules maintain their immediate set of 
neighbors for hundreds of collisional or vibrational peri- 
ods. This occurs near the temperature Ta- These local 
spatial patterns persist ever longer as the temperature is 
lowered. Interconversion between such structures occurs 
both above and below the glass transition temperature 
T g , which depends on the preparation time scale. The in- 



terconversion is called the a-relaxation when the material 
remains in equilibrium. However, when a-relaxation be- 
comes too slow and only a fraction of the interconversions 
have time to occur, the material is a glass that "ages". 
Even at cryogenic temperatures (liquid He and below), a 
certain fraction of the sample will harbor several kinet- 
ically accessible states. Interconversions can still occur 
by tunneling. These quantum motions are discussed in 
Section IV. In the final Section V, we make concluding 
remarks and highlight some open questions in the field. 

II. BASIC PHENOMENOLOGY OF THE 
STRUCTURAL GLASS TRANSITION 

Liquids exhibit a remarkable range of dynamical be- 
haviors within a relatively narrow temperature interval. 
Viscosity, for example, varies over a tremendous dynamic 
range: Fig. 2 reproduces the celebrated "Angell" plot of 
the viscosities for superooled liquids as functions of the 
inverse temperature scaled to their respective glass tran- 
sition temperatures, where the relaxation time is roughly 
one hour [6] . The temperature dependence of other struc- 
tural relaxation times, such as the inverse of the lowest 
frequency peak of the dielectric susceptibility, follow a 
similar temperature dependence and can be described by 
the so-called Vogel-Fulcher (VF) law, to a first approxi- 
mation: 



where the material coefficient D is called the liquid's 
"fragility". The Vogel-Fulcher fits work better in the 
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FIG. 3: Shown are the temperature dependences of the heat 
capacity for two substances, both for the glassy and crystalline 
counterparts. This figure is a fragment of Fig. 2 from Ref.[6]. 

vicinity of the temperature To, which is also a material 
dependent quantity. Near To, the relaxation times grow 
most rapidly and would appear to diverge, if one had the 
patience to equilibrate the liquid for cosmological times. 
Deviations from the VF law occur, of course, when the 
substance falls out of equilibrium near the laboratory T g . 

In parallel with the dynamical changes upon supercool- 
ing, there are apparent thermodynamic changes. When 
the system falls out of equilibrium, thermodynamic sus- 
ceptibilities nearly discontinuously decrease. The most 
interesting of these is the heat capacity shown in Fig. 3. 
The heat capacity allows one to monitor the entropy of 
liquid configurations of a supercooled fluid. This dis- 
continuity in the heat capacity is well approximated by 
subtracting from the measured heat capacity of the liq- 
uid the heat capacity of the corresponding crystal. By 
integrating in temperature, one obtains the part of the 
entropy that comes from the diversity of liquid configu- 
rations, see Fig. 4. This excess "configurational" entropy 
is typically of the order of a few ks per rigid molecular 
unit (such as a SiC>4 tetrahedron in silica, or an aromatic 
ring in TNB). 

When extrapolated below the glass transition temper- 
ature T g , the configurational entropy, s c , would appear 
to vanish at a temperature Tr- [9] , called the Kauzmann 
temperature. Given the magnitude of the heat capacity 
jump Ac p , at the glass transition, various ways to fit the 
temperature dependence of the configurational entropy 
can be proposed, such as [8]: 

s c = Acp(l - T K /T). (2) 

While there have been disputes about precisely how s c 
should be extrapolated [10], it is fairly clear that the dy- 
namic T) is equal to the thermodynamic Tk for all glass- 
formers, see Fig. 4. That the missing translational mo- 
tion is the main contributor to the entropy loss is explic- 
itly confirmed by the frequency dependent heat capacity 
measurements of Nagel and coworkers [11], which show 
the process contributing to the bulk calorimetry occurs 
again exactly at the time scales implied by kinetic mea- 
surements, such as the frequency dependent relaxations, 
mechanical or otherwise. That this is a time scale for mi- 
croscopic movements is also evidenced by the "plateau" 
in time-resolved neutron scattering, which confirms the 
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FIG. 4: This figure, from Ref.[8], illustrates the extrapola- 
tion of the experimentally determined configurational entropy 
(shown by squares) below the dynamic glass transition tem- 
perature T g . This graph illustrates that the Kauzmann tem- 
perature Tk is equal to the temperature To from Eq.(l), at 
which the relaxation times would strictly diverge. The circles 
represent l/ln(r), times a convenient scaling constant. 



local molecular environment rearranges on the same time 
scale [12], see also Fig. 6. 

Rather modest deviations from the consonance of time 
scales of molecular motions have excited much attention - 
these are often discussed as "decoupling." For example, 
the diffusion coefficient differs from its Stokes-Einstein 
value by about two orders of magnitude at T g . This is two 
orders of magnitude out of 14. Leaving aside these mod- 
est effects, a staggering amount of data acquired during 
the past century show the dramatic slowing is essentially 
shared by nearly all motions. 



III. CLASSICAL THEORY OF THE GLASS 
TRANSITION AND SUPERCOOLED LIQUIDS 

A. Emergence of a Free Energy Landscape of 
Aperiodic Structural States 

In gases, molecules spend relatively little time near 
each other, moving in straight lines separated by col- 
lision events, as illustrated in Fig. 5(a). In the denser, 
colder liquid state, the collision duration grows while the 
time between collisions becomes smaller. These times be- 
come comparable in the liquid state. Nevertheless, any 
two molecules are unlikely to stay together for any signifi- 
cant time. In other words, consecutive collisions near the 
critical point usually, but not always, occur between dif- 
ferent pairs, as in Fig. 5(b). The high temperature liquid 
is just a dense gas, dynamically. 
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Liquid Kinetic Regimes: 
(a) Ideal Gas 




(b) Dilute Liqaid 




(c) Landscape 




FIG. 5: The three panels demonstrate the three relatively 
distinct kinetic regimes of an equilibrium liquid, or a condi- 
tionally equilibrium liquid, as the case would be below the 
melting temperature T m . See text for explanation. The bot- 
tom portions illustrate the region of the two-particle poten- 
tials explored in the corresponding regimes. 



As temperature decreases and density increases still 
further, groups of consecutive collisions will occur ever 
more often within nearly the same set of molecules, see 
Fig. 5(c). Each molecule begins to reside within a spe- 
cific "cage" for a discernible time. Alternatively, we may 
say persistent local liquid structures form at these tem- 
peratures. Perturbative expansions, that keep track of 
the very many collision sequences within local regions, 
require computing increasingly higher order correlation 
functions [13]. At sufficiently high density, it is con- 
venient to change one's perspective of liquid state dy- 
namics from considering the dynamics of a dense gas to 
motions within sets of aperiodic crystal structures. One 
notices that even though each long-living structure in- 
volves many correlated many-particle events, each par- 
ticle within a structure is most often doing something 
quite simple. It is vibrating about a fixed location and 
only occasionally moves to a new location. This is an ex- 
perimental fact, as evidenced by neutron scattering data 
which exhibit a plateau in the time dependent structure 
factor, see Fig. 6. 

The motion of a particle is ever more confined to oc- 
cur within the cage, so individual molecular bonds dis- 
tort very little, allowing one to treat these persistent mo- 



tions as approximately harmonic at high enough densi- 
ties. When a particle takes up a new position of resi- 
dence, this implies it has moved beyond a certain thresh- 
old distance (see below) . This threshold is reminiscent of 
the Lindemann melting criterion [14] for periodic crys- 
tals. When a particle takes up a new residence position 
it usually does not do so alone. Individual bond break- 
ing and vacancy formation are just as rare processes in a 
deeply supercooled liquid as they are in crystals, which 
have comparable local density and stiffness. Instead, usu- 
ally a group of particles moves whenever a net displace- 
ment of an individual molecule occurs. 

The mostly harmonic nature of individual displace- 
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FIG. 6: The plateaus in time resolved neutron scattering 
structure function provide a direct proof of long-lived local 
structure in supercooled liquids. The curve (b) was obtained 
below the glass transition. The very flat plateau implies there 
is only a relatively small degree of structural rearrangements 
in the frozen lattice. This figure is taken from Ref.[12]: 



ments in long-living structures allows one to use a simple 
theory to describe the emergence of aperiodic crystals us- 
ing density functionals [15, 16]. In equilibrium statistical 
mechanics the free energy can be written as a functional 
of a non-uniform density. There is an entropic cost for 
forming such a non-uniform density, but if particles are 
localized, they can avoid each other to compensate. By 
combining the ideal gas localization entropy with an ef- 
fective interaction, Ramakrishnan and Yussouff [17, 18] 
wrote: 



F[p{r)} = k B T / d 3 rp(r)[l n/9 (r) 



u + i 



d 3 rd 3 r'[p(r) -p ]c(r, r'; po)W) ~ Po] + Fu 



(3) 



where F un i is the free energy of completely uniform liq- 
uid. More general functionals also have been constructed 
and used [19]. Crystallization would be described by a 
periodic p(r), but one can also use such functionals to 



examine the stability of any conceivable state with per- 
sistent aperiodic density. In view of the mostly harmonic 
nature of individual cages, we can employ the following 
variational density profile to dress any mechanically sta- 
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FIG. 7: This is a schematic graph of the free energy den- 
sity of an aperiodic lattice as a function of the effective Ein- 
stein oscillator force constant a (a is also an inverse square 
of the localization length used as input in the density func- 
tional of the liquid. Specifically, the cureves shown character- 
ize the system near the dynamical transition at Ta, when a 
secondary, metastable minimum in F(a) begins to appear as 
the temperature is lowered. 



ble zero temperature structure: 

p(r) ee p(r,{r t }) = £ e -(--<) 2 . (4 ) 

i 

Here {r^} represent the atomic coordinates of a particu- 
lar, but generic aperiodic lattice characterized by some 
average density n = 1/a 3 , where a is the average lat- 
tice spacing. The variational parameter a is a conve- 
nient measure of localization within such a many-particle 
cage. While nonzero a characterizes a localized regime, 
where transient local quasiharmonic environment forms, 
the same function is also capable of describing the com- 
pletely dclocalized regime of the uniform liquid by tak- 
ing a = 0. Singh, Stoessel, and Wolynes [16] estab- 
lished that at sufficiently high density, at a given tem- 
perature called Ta, the free energy from Eq.(3) devel- 
ops a metastable minimum as a function of a, see Fig. 7. 
The transition is very similar to a spinodal crystalliza- 
tion, but to an aperiodic structure. Similar results were 
obtained earlier by Stoessel and Wolynes [15] using a self- 
consistent phonon approach. The density functional and 
self-consistent phonon approaches both suggest that at 
Ta, metastable structures form. This metastability indi- 
cates there would be a corresponding barrier demarcating 
one metastable minimum in the entire set of such min- 
ima from the uniform liquid state. The emergence of 
the minimum appears as a first order transition, where a 
plays the role of an order parameter, whose value changes 
discontinuously, during the transition, from a = to 
a = ao- 

The quantity measures the extent of vibrational 

motion within a local metastable structure. Motions of 
any greater extent, at the mean-field level, signify a tran- 
sition between thermodynamically distinct states. The 
quantity 1/^/ao is the vibrational amplitude at the me- 
chanical stability edge, fully analogous to the famous Lin- 
demann length for ordinary crystals [14]. We will denote 



it as dh'- 

d L = (5) 

The Lindemann ratio /a is roughly one-tenth for all 
crystals [20-22], and a similar value emerges from the 
microscopic calculations for aperiodic structures. This 
value is also consistent with the magnitude of the plateau 
in the neutron scattering time correlation function [12]. 
While in terms of ao the transition is first order, there 
are many aperiodic states randomly distributed in free 
energy. We therefore call this transition the "Random 
First Order Transition" (RFOT) because, while first or- 
der in a, it results in forming many random (infinite life- 
time) "phases" which arc distinct both morphologically 
and spatially. Because of the multiplicity of states there 
is no latent heat but instead at the mean field level, there 
would be a heat capacity discontinuity much as occurs 
when a liquid falls out of equilibrium. 

Even though the metastable minimum of F(a) is 
higher in free energy than the uniform liquid state, we 
realize that the F(a) curve was computed for a single 
(typical) aperiodic lattice. The full liquid free energy 
available at a incorporates the multiplicity of all such 
(typical) aperiodic states, i.e. e ScN for a region encom- 
passing N molecules. The multiplicity of aperiodic states 
at the "spinodal" temperature Ta obviously must imply 
the localized minimum is metastable. Whether at a suf- 
ficiently low temperature it can reach F un i and thus im- 
ply the existence of an observable thermodynamic state, 
is not immediately obvious. Yet several rigorous argu- 
ments suggest this is indeed the case. Kirkpatrick and 
Wolynes [23] showed that the mode-coupling viscosity 
catastrophe, which can be obtained by summing rccolli- 
sion events, also corresponds to the same transition at Ta 
predicted by self-consistent phonon and variational ap- 
proaches. They noticed the connection with Potts spin 
glasses. Kirkpatrick and Thirumalai clarified this rela- 
tion which KW had explicitly showed only in high dimen- 
sion for fluids [24, 25]. Like the approximate treatment 
of the fluid system, the exactly soluble mean-field Potts 
glass also exhibits a dynamic transition at a tempera- 
ture Ta above its thermodynamic glass transition. Kirk- 
patrick and Wolynes [26] went on to show that these infi- 
nite range Potts spin models, too, exhibited a Kauzmann- 
like entropy crisis at the lower thermodynamic tempera- 
ture Tk- They also suggested that for finite range Potts 
spin glasses and supercooled liquids the Ta transition 
would be smeared by droplet-like excitations driven by 
the configurational entropy. These "entropic droplet" ex- 
citations would provide the route to equilibration below 
the mean field dynamical temperature that corresponds 
to the mode coupling transition. Their argument also ex- 
plained why the dynamic To and Tk should be the same. 
These entropic droplet excitations also form the basis of 
the microscopic developments of RFOT theory that in re- 
cent years, have explained supercooled liquid and glasses 
quantitatively. We will describe entropic droplets in a 
somewhat different language in the following section. 
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The said mean field picture, which employed density 
functional ideas to treat the fluid aspects and the analogy 
to the exactly solvable Potts spin glasses, has been de- 
veloped more formally using the methods of replica field 
theory. By making reference to the molecular fluid con- 
figuration at one time as a fiducial structure, the replica 
method can be used to construct stability criteria and 
develop approximations for the configurational entropy 
directly [27]. Numerous analyses using these tools con- 
firm the earlier mean field developments. These more 
explicit methods have also been tested near the dynamic 
transition by computer simulations [28]. The reformu- 
lation of the theory using replicas also has allowed the 
activated events to be studied in a systematic way when 
the range of the interactions is finite but large [29-31]. 
These mathematically controlled developments increase 
the confidence in the simpler constructive arguments for 
structural glasses that we will describe below. Another 
direction in which the microscopic calculations can be 
taken is to explore the molecular origins of the con- 
figurational entropy. Hall and Wolynes recently used 
these approaches to explain why network materials with 
more constraints (like SiC^) behave like "strong" liquids 
while simple van der Waals systems without bonding con- 
straints are more fragile [32] . 

It is important to note that within the RFOT theory, 
detailed microscopic calculations starting from the inter- 
molecular forces can in fact be carried out. Presently, 
only the mode coupling theory can make a comparable 
claim. Mode coupling theory (MCT) is not an orthogonal 
approach, however. As Kirkpatrick and Wolynes pointed 
out in 1987 [23] the two approaches are effectively equiv- 
alent near the dynamical crossover. Recent successes of 
MCT in describing re-entrant glass transitions in attrac- 
tive colloids, etc. thus buttress the general RFOT picture 
from a microscopic viewpoint. Microscopic calculations 
set the stage for the RFOT theory. Much as in think- 
ing about conventional phase transistion, it is important, 
however, to distinguish the main RFOT ideas from the 
liquid state engineering details for specific systems. It is 
clear, for example, that the deep supercooled regime of 
van der Waals liquids probes intermolecular forces at a 
more intimate length scale than near the critical point. 
Time-honored prescriptions based on perfectly hard po- 
tentials thus will break down. This breakdown is most 
clearly evidenced by the experimental fact that the iso- 
choric activation energy of the viscosity is not zero be- 
low Ta, while it is nearly so above that point. Just as 
for studies of crystals [33], the popular Weeks-Chandlcr- 
Andersen prescription [34] , for dividing the potential, re- 
quires re-examination in deeply supercooled liquids. 



B. The Library of Local States 

In mean field theory, the lifetime of an aperiodic struc- 
ture is infinite. But in reality the lifetime of a system 
with finite range interactions is finite because rearrange- 



ments can occur independently within regions of finite 
spatial extent. The number of relaxation events per unit 
time scales with the system size. A mole of water will 
undergo roughly 10 39 barrier crossing events a second. 
These events arc local in character in that they affect 
each other little beyond some critical distance. 

The locality of structural processes in supercooled liq- 
uids is intrinsic in the entropic droplet concept. Recently 
their quantitative aspects have been constructively es- 
tablished within the RFOT theory [35, 36]. Our present 
discussion will be based on the "library" construction 
of local configurations [37] , which makes the ideas more 
transparent than in earlier discussions and also allows 
the treatment of the "aging" regime. This library con- 
struction is essential for the latter far-from-equilibrium 
case. An elegant formal analysis along similar lines has 
been made by Bouchaud and Biroli [31]. This construc- 
tion can also be quantized to deal with the phenomena of 
cryogenic temperatures [38]. In the library construction, 
one first averages over the vibrational modes of the super- 
cooled liquid. These equilibrate generically in a time less 
than a picosecond. The variational ansatz used in den- 
sity functional theory from Eq.(3) operationally defines 
such an averaging. The lifetimes of the resulting aperi- 
odic crystal states depend on the size of the region. For 
the library construction to be strictly valid, the processes 
must be slower than the sound modes of similar length 
by three orders of magnitude or so. Fortunately, this cri- 
terion covers most of the dynamical range accessed by 
supercooled melts, except perhaps for the high frequency 
processes referred to as the Boson peak in liquids [39]. 
Transitions between the aperiodic crystal states defined 
above give rise to the configurational entropy. 

To convert from one state to another state we must 
consider not only direct paths connecting one structure 
to a neighboring one; but most importantly, one must ex- 
amine all thermally realizable dynamically connected se- 
quences containing many structural states. By "dynami- 
cally connected" we mean the trajectory is a sequence of 
processes with small barriers that involves a translation 
of a single bead by a distance not exceeding the Linde- 
mann length cLl- The lowest energy path leading to a 
thermally representative liquid state will define the most 
probably escape trajectory from the initial state. The 
highest free energy point along the path will determine 
the bottleneck, or the critical value of the progress coor- 
dinate. The energy at the bottleneck will determine the 
barrier height and will thus give the activation part of 
the escape rate. A very nice example of a computational 
study of such a sequence is provided by Saksaengwijit 
and Hcuer [40] . They show that even in high temperature 
silica, whose rates are nearly Arrhenius, reconfiguration 
events proceed through several steps, see Fig. 9. 

Now, the totality of the aperiodic states that could be 
constructed with a region of fixed boundaries forms a "li- 
brary" of local states. This construction is illustrated in 
Fig. 8. In this figure, the circles do not signify the instan- 
taneous positions of the molecules but rather their most 



7 



o o o u 
o ooo° u H 

oo ooo°o 
O n u o O ° Q 

Initial Configuration 
Global Library 



o o o u 
o ooo° u H 

rPooo °c 

O n u n O O 

N=S 

Local Library N=5 



000 u 
o ooo° u H 

On u n O °0 

Local Library N=7 




FIG. 8: This figure is taken from [37]. In the upper panel on 
the left a global configuration is shown, chosen out of a global 
energy landscape. A region of N = 5 particles in this config- 
uration is rearranged in the center illustration; this subset is 
taken from another aperiodic state or from a different loca- 
tion of the same liquid state. The original particle positions 
are indicated with dashed lines. A larger rearranged region 
involving N = 7 particles is connected dynamically to these 
states and is shown on the right. In the lower panel, the left 
most figure shows the huge density of states that is possible 
initially. The density of states found in the local library origi- 
nating from a given initial state with 5 particles being allowed 
to move locally is shown in the second diagram. These ener- 
gies are generally higher than the original state owing to the 
mismatch between the two structures. The larger density of 
states where 7 particles are allowed to move is shown in the 
right most part of this panel. As the library grows in size, 
the states as a whole are still found at higher energies but 
the width of the distribution grows. Eventually with growing 
TV, a state within thermal reach of the initial state will be 
found. At this value of N* we expect a region to be able to 
equilibrate. 
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FIG. 9: Even in "strong" liqiuds, reconfigurations occur 
through a series of steps as in the library constuction. Sak- 
saengwijit and Heuer have delineated these events in silica. 
Their results are shown with their Fig. 4 [40] . Panel (a) shows 
the temperature dependence of the average waiting time in 
the low energy inherent structures of simulated silica. In (b), 
a specific escape trajectory from an inherent structure is rep- 
resented. The percentages indicate how often such a structure 
recurs to its initial state. Panel (c) gives the distribution of 
escape barriers from IS. 



labelled by a dummy index j whose bulk energy is: 



bulk 



Ej — TS V 



br- 



(7) 



Choose a boundary B encompassing N beads some- 
where within the initial structure, and remove its con- 
tents. We then cut an identically shaped region out of 
the structure j and paste this region in the void in the 
initial structure. Call the modified region of the initial 
state a "droplet" . The free energy of this trial configura- 
tion is: 



bulk 



bulk 



(8) 



probably positions after vibrational averaging. Further, 
the circles correspond to effective individual structural 
units, or "beads" . Chemical intuition often makes clear 
how to decompose conceptually a molecule into beads. 
A rigid chemical group, a single benzene ring, or a com- 
pact side-chain in a polymer will usually constitute an 
independently moving unit. The volumetric bead den- 
sity can be quantified unambiguously by comparing with 
the fusion entropy, if the corresponding crystal exists, 
since bead motions are also frozen in the crystal; see the 
detailed discussion in [41]. At any rate, a bead is usually 
a few angstroms in size. 

Consider an initial configuration spanning some ex- 
tended region. Its "intrinsic" bulk energy can be written: 



fl^bulk 771 



-TS- 



vibri 



(0) 



which combines the energy proper of the particular ar- 
rangement of the most probable molecular coordinates 
with the free energy of the high frequency vibrations. 

This state may move to some other thermally repre- 
sentative aperiodic crystal state over the same region, 



This free energy includes a mismatch penalty Vj^ n be- 
cause no matter how carefully we placed the surface B in 
the structure j, the contents of the latter region were vi- 
brationally averaged with a constrained surrounding dif- 
ferent from the one in the initial state and thus some 
"elastic energy" must be paid. The alternative state 
inside the droplet may be thought of as the original 
molecules somewhat displaced from their original posi- 
tions. The cut-and-paste routine realizes a finite number 
of thermodynamically distinct droplet states. This fi- 
nite number is given by the configurational entropy, and 



scales exponentially with the droplet size: e 



s c N 



In the 



central panel of Fig. 8, we show, for illustration, just the 
subset of droplet states available at size N = 5. This 
subset would constitute an astronomically minute frac- 
tion of the whole set of states available to an extended 
sample. Now, consider the droplet states at different, but 
close sizes. Some states at both sizes will be structurally 
very close, implying high degree of dynamical connec- 
tivity. Sequences of such dynamically connected droplet 
states constitute escape trajectories from a given initial 
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extended state to some other extended aperiodic crystal 
state. Because one may associate with any of those se- 
quences of droplet states a sequence of droplet sizes, one 
may therefore speak of nucleation of one aperiodic crys- 
tal state within another and of the corresponding prop- 
agation of a domain wall separating those two aperiodic 
crystal states. In this way, the mismatch penalty may be 
thought of as the domain wall tension. 

Because of the surface penalty, the lowest energy, per 
particle, in a droplet library is likely to be well above 
the energy, per particle, in the initial extended state. 
Still, the energies within such finite subsets are always 
distributed, roughly in a Gaussian fashion. The lowest 
energy in the subset, it turns out, is determined by a com- 
petition between two factors: On the one hand, with the 
growing droplet size, the mismatch penalty will increase 
thus shifting the bulk of the distribution of the subset 
energies upwards. On the other hand, the weight of the 
distribution, i.e. e s " N , increases exponentially with N, so 
that the lowest energy state will be found at an increas- 
ingly lower position. Eventually, at large enough droplet 
size N*, the latter trend will compensate the former, so 
that the droplet will be guaranteed to have the energy 
per particle equal to that of an arbitrarily extended state. 
Thus to a region of size N* or larger, all thermodynami- 
cally relevant bulk liquid states will be available; i.e. such 
a region is expected to be typical of equilibrium at the 
temperature T. A supercooled liquid is a mosaic of aperi- 
odic crystals! In the following we will put this qualitative 
discussion in formal, quantitative terms. 

Mismatch Penalty between Aperiodic Crystal 
States. To a first approximation from the density func- 
tional viewpoint, computing the mismatch energy is sim- 
ply a matter of counting up the missing interactions from 
some unsatisfied local contacts in an interface region. 
The mismatch penalty would thus scale with the inter- 
face area itself, times the energy of the unsatisfied bond: 
T cx er d_1 Z m i cro , where e and r are the interface energy 
density and the region size respectively [26] . This is also 
the result of instanton calculations based on replica meth- 
ods [42, 43]. However, the situation is a bit more com- 
plicated as we must recognize that some of the states 
in the library will in fact match much better than do 
others because they are partially random. The scaling 
of thermal averaged T with TV will thus be weaker than 
expected for interfaces between entirely distinct phases. 
The fraction of the better matching configurations is sig- 
nificant enough to actually partially short-circuit the con- 
ventional surface tension, renormalizing it to a smaller 
value. From an elastic theory point of view, we expect 
the displacement fields due to the local structural rear- 
rangements in our aperiodic structures not to exhibit long 
range correlations. Accordingly, the deviation of total 
strain within a region of size M from its average value 
would scale with the usual y/M, and is of cither sign 
with equal probability. This means one may generically 
impose a random external field with the usual Gaussian 
statistics to lower the energy of a region of size M by an 



amount scaling with \[M ' . A curved interface between 
two states will distort so as to lower the local free energy 
to take advantage of this random contribution. The en- 
ergy compensation will scale, again, as the square root of 
the variation of the volume occupied by either of the two 
phases, due to the boundary distortion. The final shape 
of the interface will be determined by a competition be- 
tween this stochastic energy compensation and the cost 
of increasing the area of a flat interface. Consistent with 
this, the scaling of surface tension with r differs from the 
conventional r^ d ~^/ d . This situation is analogous to the 
problem of the interface between the spin-up and spin- 
down regions in the Random Field Ising Model (RFIM), 
which has been treated by Villain [44]. This analogy 
was explicitly exploited by Kirkpatrick, Thirumalai, and 
Wolynes (KTW) [35] in deducing how the droplet inter- 
face tension scales with the droplet size. (One may in 
fact show that the two problems can be mapped onto 
each other, with the error not exceeding the discrepancy 
between the microcanonical and canonical averages [45].) 
A pleasant dividend of this argument is that the hyper- 
scaling relation connecting the heat capacity and the mo- 
saic length scaling near Tk ■ (2 — a — vol), - is restored 
[35]. The KTW argument suggests the effective surface 
tension coefficient is renormalized according to: 

a(r) ~ aoir/a)-^- 2 ^ 2 , (9) 

where do is the surface tension coefficient at the molec- 
ular length scale. In all spatial dimensions d > 2, the 
scaling of the resulting mismatch free energy 

T oc r d ' 2 (10) 

is numerically inferior to the conventional r( d_1 )/ d scal- 
ing, but of course this is quantitative only for sufficiently 
large interfaces. It also follows that the thickness of the 
distorted interface scales with the radius itself [35, 44] 
consistent with our a priori expectation that the mis- 
match energy is determined by the subsets of static struc- 
tures that smoothly interpolate between two generic ex- 
tended aperiodic states. It may be said that an inter- 
face between two thermally averaged aperiodic crystals 
is, strictly speaking, not thin but is always "wetted" by 
other states. 

On short length scales, the coefficient a in Eq.(9) must 
correspond to the mismatch penalty at the molecular 
length scale. This suggests a way in which tro can be 
estimated in zeroth order, at the molecular length-scale. 
The mismatch requires the particle at the interface to still 
be within the Einstein oscillator localization volume d\, 
instead of its volumetrically available space a 3 . Yet in the 
free energy functional it does not receive the full benefit 
of its neighbors' being localized (only half of them arc) 
and thus staying out of the way. Since these free ener- 
gies must balance at Tk, there is an additional mismatch 
penalty on this scale that is related to the localization en- 
tropy cr ~ ±k B T\n(a 3 /d 3 L ) = (3/4)fc B Tln(a 2 /d|). This 
expression for ctq contains all the essential parameter seal- 
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ing, but obviously there may be some numerical uncer- 
tainty in applying the long distance scaling of cr(r) all the 
way down to the shortest lengths. Nevertheless, by virtue 
of the slow, logarithmic dependence on the constants, 
along with the near universal value of the Lindemann ra- 
tio di/a = 0.1, this argument suggests Co is a universal 
multiple of the glass transition temperature itself. This 
is the key to the later quantitative results obtained in 
the RFOT theory. Using the density functional approxi- 
mation of Eqs.(3) and (4), Xia and Wolynes specifically 
obtained the surface tension at the molecular scale [36]: 



(JO = T- 



3k B T 



ln[(a/rf L ) 2 Ae], in a(r) = cj^a/r) 1 ' 2 . (11) 



This value works very well quantitatively for real ma- 
terials. The chemical universality of <To reminds one of 
Turnbull's rule used in treating crystallization that em- 
pirically states the tension between a periodic crystal and 
its melt is a universal multiple of T m [46]. It is conve- 
nient to write the mismatch free energy in terms of the 
number of reconfigured particles: N = (47r/3) (r/a) 3 : 



T = jVN, where 



2V3tt 



7 



k B T\n 



{a/d L f 



7re 



• (12) 



While the mismatch energy T in j is doubtless distributed, 
as long as we are dealing with typical liquid states, we 
can use its typical value T. 



log of the multiplicity, S c (^g lk ), is nothing but the con- 
figurational entropy at equilibrium, which can be mea- 
sured by calorimetry: ^(^j^ 111 *) = S c . 

The N dependence of all the parameters in Eq.(13) 
becomess obvious from Eq.(8) upon recalling that S c 
and <I> bulk are bulk parameters and scale with N itself: 
S c = Ns c , $ bulk = N(/) hulk . One thus obtains a sim- 
ple expression for the typical free energy profile during 
a structural rearrangement after optimizing with respect 
to the saddle point energy. It varies with size as: 



F(N) = [0 e b 9 ulk (T) - ^ k (T)]N + 7 ViV - Ts c N. (14) 

Here, we have used Eq.(8). F does reach down to arbi- 
trarily low energy states, but to escape, one needs to pass 
first over the maximum, which gives the typical barrier 
F$ for reconfiguration from a given initial state. When 
F = 0, a typical state has already been reached. 

Consider supercooled liquids equilibrated above the 
glass transition temperature T g . The initial liquid state 
is thermodynamically typical of the temperature T thus 



ibulk 



(T) = (f>^ lk (T). The resulting nucleation profile is 



quite simple: 



F(N)\ 



T>T g 



= jy/N - Ts c N. 



(15) 



Clearly, structural transitions are driven by configura- 
tional entropy alone! Eq.(15) immediately gives an in- 
verse scaling of the most probable relaxation barrier with 
the configurational entropy density that automatically 
yields the Vogel-Fulcher law: 



C. Activated Motions between Local States 

We see, via the library construction, that if too small 
a region is reconfigured, even the lower energy paths will 
have a monotonically increasing energy with the size. 
But as the region is made larger, some fraction of paths 
will curve down in energy. At sufficiently large size, there 
will always be a trajectory ending at an energy within the 
thermally relevant liquid energy range. In computing the 
typical escape rates from a liquid state, one therefore only 
needs to escape by passing through droplet configurations 
with the size corresponding to the lowest flux. Formally 
averaging the escape flux over the ensemble of transition 
state droplets yields [37]: 



k 



= r 1 jWi/cde 8 ' 



{^)/k Bp -{4>r 



4> in )/k B T 



r -l e Sc(*™t )/fcB e -Wopt-0/fcBT 



= '0 



r i e Sc(*^ lk )/fc Be -(^ 9 -^; b )/fc B T 



(13) 



The intrinsic activation free energy represents a critical 
droplet-configuration energy from Eq.(8) and the factor 
e Sc(* :j )/k B gjvgg t ne multiplicity of liquid configura- 
tions at the energy <f>j. (c^ is a normalization constant.) 
The energy <f> op t, maximizing the integral, must be as- 
signed the equilibrium energy value. Consequently, the 
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As c T AAc p (T - Tj 



K j 



(16) 



Here we used the specific form (2) for the configura- 
tional entropy. Note that if it were not for the surface 
tension rcnormalization, see Eq.(10), the s c dependence 
in the denominator of the middle expression would be 
quadratic, not linear. The inverse scaling of the relax- 
ation barrier with the configurational entropy was de- 
rived by KTW [35] , but was proposed originally by Adam 
and Gibbs (AG) [47], who stipulated that there be a 
smallest rearranging unit in a liquid. This AG unit would 
be characterized by two configurations, whose size was 
assumed independent of temperature. The Adam-Gibbs 
argument does not reconcile how the existence of such a 
special fixed size is compatible with that argument's use 
of a temperature dependent, extensive configurational en- 
tropy. In contrast, within RFOT theory the length scale 
of activation is determined by the underlying Hamilto- 
nian and varies with temperature. The RFOT thoery 
precisely predicts the way in which the critical size of the 
nucleation barrier scales in Ref. [35] : 

(17) 



oc 



{Nt)1/3 a * {T-T K )W 



This scaling law is consistent with a specific heat dis- 
continuity and the usual hyperscaling relation for contin- 
uous transitions [35] . RFOT theory directly shows the ki- 
netic and thermodynamic anomalies in supercooled melts 
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are intrinsically related. Furthermore, given the value of 
(To from Eq.(ll), Xia and Wolyncs established an amaz- 
ingly simple relation between the kinetic fragility and the 
heat capacity discontinuity: 

D = 32./Ac p . (18) 

The universal numerical constant 32. comes from the nu- 
merical constants in the argument and the universal value 
of the Lindemann ratio. Ac p is the heat capacity jump 
per bead. Owing to the universality of the Lindemann ra- 
tio, RFOT theory implies the material dependence of the 
fragility comes essentially from the heat capacity only. 
Since the Lindemann length enters under the logarithm, 
its small variations do not affect the constant signifi- 
cantly. In any event, a recent argument by Lubchcnko 
[48] shows the Lindemann ratio di/a would be expected 
to vary at most my 10% in supercooled liquids. 

The present argument assumes the liquid to be suffi- 
ciently deeply supercooled so that deviations from the 
activated picture are small. This is true already when 
relaxations are about three orders of magnitude slower 
than the times of vibrational temperature equilibrium. 
The activated picture thus applies at t > 10~ 9 sec, 
which covers most of the liquid dynamic range. At 
higher temperatures, collisional effects on the viscosity 
become important. The agreement of the relation of 
Eq.(18) between D and Ac p with available data is im- 
pressive. Nevertheless, to avoid fitting ambiguities, it 
appears best to use the measured fragility index near T g : 
m = T[dlog 10 T(T)/d(l/T)]. The index m scales roughly 
inversely proportionally with D, but is less dependent on 
the fitting ambiguities arising from the crossover to the 
high temperature region. Like D, m also follows from 
Ac p without adjustable parameters. Fig. 6 shows how 
theory and experiment compare for the predicted m from 
thermodynamics and the measured kinetic values. There 
are a small number of outliers. In examining those excep- 
tions one must bear in mind, however, the RFOT theory 
applies strictly only to purely amorphous and fully equi- 
librated melts. Deviations from the RFOT predictions 
are expected when samples exhibit partial crystallinity, 
often present in polymers, or other types of local order, 
as in decalin or in alkali borates (see e.g. [49]), 

Whether such local order effects exist may actually be 
judged, to some extent, using another simple result that 
follows from the RFOT theory: Near T g , relaxations are 
strictly activated, implying the typical relaxation bar- 
rier relative to the temperature depends only logarithmi- 
cally on the relaxation rate: (F$ /T)t 9 = ln(r/ro). The 
quenching rate in the laboratory is limited by the sam- 
ple's heat conductance, on the faster side, and the exper- 
imenter's patience and stable temperature conditions, on 
the longer side. Combining a one hour time scale with 
Eqs.(12) and (16), RFOT predicts the configurational en- 
tropy at T g , per bead, should be [36]: 

s c (T g ) ~ 0.8, (19) 




m, theory 

FIG. 10: The horizontal axis shows the value of fragility as 
computed from the thermodynamics by the RFOT theory, 
and the vertical axis contains the fragility directly measured 
in kinetics experiments. Here m is the so called fragility in- 
dex, defined according to m = T[dlog 10 r(T)/d(l/T)]. m is 
somewhat more useful than the fragility D, because devia- 
tions from the strict Vogel-Fulcher law, r = T e DTK ^ t ^ Tk \ 
are often observed, see text, m essentially gives the appar- 
ent activation energy of relaxations at T g , in units of T g , it 
is roughly an inverse of D. In evaluating m theoretically, one 
needs to know the size of the moving unit, or "bead", in each 
particular liquid. The latter can be estimated using the en- 
tropy loss at crystallization and scaling it to its Lennard- Jones 
value [41], resulting in m thcor oc A g ^f(T°) x ^ah^ 9 > in 
view of the near universality of s c (T g ) (see text). This figure 
is taken from [50]. 

latter time ratio, at T g , is very large: 10 16 or so, ren- 
dering the estimate above nearly universal. A brief look 
at the configurational entropies of polymers, such as in 
Fig. 4 of Roland at el. [51] immediately gives away there 
are many partially crystalline polymers. For example, 
PVC has s c (T g ) ~ .1 clearly indicating a large degree of 
crystallinity, consistent with X-ray determinations of the 
latter. 

Another important universality predicted by the 
RFOT theory is the size of the cooperatively rearrang- 
ing region at the glass transition temperature. For the 
droplet size N > N*, where F(N*) = 0, one expects all 
thermally relevant liquid states have become available; 
there is no further "growth" of an alternate phase. Typ- 
ically then TV* beads are reconfigured during a typical 
relaxation event. Using 

N* = (£/a) 3 , (20) 

one easily finds [36]: 



depending only logarithmically on the ratio t/tq. The 



£| t =t s = 5.8 a, 



(21) 
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for t/to = 10 17 , i.e. at the relaxation time of the order an 
hour. This coopcrativity length scale £ is not structural, 
but dynamic, and hence must be probed by nonlinear 
dynamic experiments. Because any region of size £, but 
no smaller, may relax, we can regard a supercooled liquid 
as a "mosaic" [36] of cooperatively rearranging regions 
or "entropic droplets" . Like the critical radius , the 
dynamical heterogeneity of length scale £ scales according 
to 



(T-T K )V^ 



(22) 



Eqs.(21) and (22) are arguably the most specific micro- 
scopic predictions of the RFOT theory; they are certainly 
consistent with the earlier semi-quantitative results of the 
4D exchange NMR experiments by Tracht at el. in PVAc 
[52]. Concurrently with the publication of the theory or 
a bit later, other types of non-linear studies have been 
performed on several supercooled liquids which all quan- 
titatively confirm the RFOT prediction. These studies 
include experiments using nanometer dielectric probes by 
Russell and Israeloff [53] . Deviations of the hydrodynam- 
ics of small probes from the Stokes-Einstein relation [54] 
are also consistent with this length scale as earlier pre- 
dicted by Xia and Wolynes [55], who showed deviations 
from the Stokes-Einstein relation are expected for probes 
smaller the length £. These deviations are especially no- 
table in that they explicitly show the cooperativity length 
is temperature dependent [55]. Very recently, evidence 
of a dynamic cooperative length, which increases with 
lowering the temperature, have been obtained using a 
rigorous inequality based on the non-linear susceptibil- 
ity. Berthier at el. [56], again, decisively confirm the 
result in Eq.(21). Finally, while the detailed tempera- 
ture dependence of the cooperativity length is somewhat 
modified by the mentioned barrier softening effects, the 
universality from Eq.(21) is robust, in the whole range of 
liquid fragilities [41]. 

Eq.(15) only gives the typical droplet nucleation pro- 
file. Specific activation paths will reflect local variations 
of the liquid landscape, i.e. they depend on the initial 
local configuration. These variations are encoded in the 
local density of states (DOS) that can be connected to 
the original configuration. A higher than average density 
of states implies more configurations are available, lead- 
ing to a lower activation free energy barrier and a smaller 
number of molecules participating in a structural transi- 
tion, see Eqs.(16) and (17). Local variations in the DOS 
correspond to fluctuations of the configurational entropy. 
Entropy fluctuations, on the other hand, follow from the 
standard formula ((A5 C ) 2 ) = ACp. Only the configu- 
rational part of the heat capacity, as given by the heat 
capacity jump at the glass transition, enters. Just this 
logic was followed in Rcf . [57] . At this level of approxima- 
tion, the entropy fluctuations, which arc approximately 
gaussian (surely, for N* > 6), lead to a nearly gaussian 
distribution of barriers, via Eq.(16). The ratio of the 
barrier distribution width to its most probable value de- 
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FIG. 11: This figure shows the correlation between the liq- 
uid's fragility and the exponent (5 of the stretched exponential 
relaxations, as predicted by the RFOT theory, superimposed 
on the measured values in many liquids taken from the com- 
pilation in Ref. [58]. The dashed line is based on a simple 
gaussian barrier distribution, with the width mentioned in the 
text. The solid line takes into account the averaging effect of 
environmental rearrangments surrounding a mosaic cell, so 
that the barrier distribution to the right of the most probable 
value is replaced by a narrow peak of the same area; the peak 
is located at that most probable value. This figure is taken 
from [57]. 



pends only on the liquid's fragility [57] : 
SF* 1 



mp 



(23) 



As a result of the barrier distribution, bulk measurement 
of relaxations will produce non-exponential time decay, 
or, non-Lorentzian profiles in the frequency domain. The 
former is often fitted with a stretched exponential: 



p(t) = e 



-(t/toV 



(24) 



Assuming the barrier distribution is purely gaussian, one 
indeed recovers a decay profile closely resembling the 
stretched exponential form, where the corresponding ex- 
ponent P is related to the fragility in the following simple 
manner [57]: 



/? ~ [l + (5F*/k B T) 2 ] 



'1/2 



(25) 



This correlation is shown by the dashed curve in 
Fig.ll. While already qualitatively consistent with ex- 
perimental data, further improvement in agreement is 
immediately achieved upon realizing that if a particu- 
lar region happens to be relatively short on alternative 
structural states, a slightly distinct, and hence overlap- 
ping (!), region is not likely to be so disabled. Since the 
environment of the original "slow" domain will change 
such a slow region will relax sooner than expected. Xia 
and Wolynes suggested therefore the barrier distribution, 
for a-relaxation, should be asymmetric, with relatively 
less weight on the high barrier flank. The simplest way 
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to incorporate the above considerations into the theory 
is to use the original gaussian distribution, but with the 
slow half of that distribution replaced by a delta func- 
tion carrying the same total weight, and centered at the 
most probable value of the original gaussian. This ap- 
proximation introduces no new parameters. Using this 
simple form still leads to a correlation of (3 with D and 
reproduces the empirically known correlation between (3 
and D quantitatively [57], see the solid line in Fig. 11. 
The RFOT theory clearly predicts the degree of non- 
exponentiality is temperature dependent, consistent with 
findings of Dixon and Nagel [11]. 



D. Dynamics near the Crossover 

The temperature Ta is a mean-field spinodal where 
one expects a siginificant barrier softening, as Ta is ap- 
proached from below. Above Ta transport is no longer 
strictly activated. Instead, weakly correlated collisions 
become the dominant contributor to the liquid's viscous 
response. The barrier softening effects have been quanti- 
tatively assessed by us [41]. These effects are particularly 
important at viscosities below 10 Poise or so, for all con- 
sidered substances. 

For salol, the theoretically derived contribution to the 
relaxation stemming from activated processes is shown 
as the solid line in the left pane of Fig. 12. The theory 
uses experimentally derived values of the configurational 
entropy and the high temperature limit of the relax- 
ation times using the prefactor given by collision theory. 
The softening correction computed by Lubchcnko and 
Wolynes (LW) relies on the idea that at small droplet 
radii, there is less interface wetting. The correction 
terms in their analysis involve only one fitting param- 
eter, which is the spinodal dynamical crossover temper- 
ature Ta, where the mean- field barrier between liquid 
states would vanish. Clearly, at viscosities less than 10 
Poise or so, where the theoretical and the experimental 
curve bifurcate, the transport becomes largely collisional. 
For several systems, the temperature where the experi- 
ment and theory diverge, turns out to lie rather close to 
the temperature T c , where Stickel at el. [59] have found 
a "kink" in the temperature dependences of relaxation 
times. Often, two distinct VF forms have been used to 
fit the data in the two temperature ranges, separated by 
T c [59]. 

While the LW analysis takes the temperature Ta from 
fits, a recent argument of Stevenson, Schmalian, and 
Wolynes [60] provides a microscopic description of the 
onset of activation-less reconfigurational motions. The 
overall high free energy profile, Eqs.(14) and (15), comes 
from the relative scarcity of phase space trajectories that 
are sequences of low cost local moves. The overall free 
energy cost should not only include competition between 
the mismatch energy of mobilizing particles against im- 
movable neighbors and the entropy gain so achieved, but 
also the number of different ways a connected pattern 
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FIG. 12: Experimental data (symbols) for salol's viscosity 
[59], superimposed on the results of the fitting procedure 
(line) from [41] are shown. Ta is the temperature at which 
the meanfield barrier vanishes, indicated by a tickmark. The 
configurational entropy, used in the fitting procedure, was 
extracted from experiment, see Ref.[8] and Eq.(2). The tem- 
perature T cr signifies a cross-over from activated to collisional 
viscosity, dominant at the lower and higher temperatures re- 
spectively. The temperature is varied between the boiling 
point and the glass transition. The r.h.s. pane depicts the 
temperature dependence of the length scales of cooperative 
motions in the liquid. The thick solid and dashed lines are r* 
and £ respectively. This figure is taken from [41]. 



of mobile particles can be placed on an aperiodic lat- 
tice. Keeping track of such patterns can be done, for in- 
stance, by counting contiguous percolation clusters [61] 
or by counting impenetratable strings emanating from a 
common origin. Here, again, the probability to find a 
contiguous cluster grows with increasing configurational 
entropy The estimates based on percolation clusters or 
strings give distinct but comparable values of the critical 
configurational entropy at which non-activated reconfig- 
uring short-circuits become possible [60]: s^" = 1.28fc_B 
and s^ tring = 1.13/cs. The corresponding temperature at 
which a crossover to non-activated transport occurs, is 
given by 



T K 



T 



K 



Ac p Tg 



(26) 



A similar result is obtained for the "string" transition. 
Here, the functional form for the configurational entropy 
from Eq.(2) has been used. In Fig. 13, we show the the- 
oretically derived values of T c pGrc and T^ tring , in compar- 
ison with the Stickcl's T c . Note that no adjustable pa- 
rameters are used in Fig. 13. It is interesting that the 
mathematics of the string calculation is isomorphic to 
the Hagedorn transition in the string theories of particle 
physics [62]. 

We see that RFOT theory predicts deviations from the 
strict Vogel-Fulcher dependence of the relaxation times 
on the temperature. These deviations are expected to 
occur also at the faster side of the full dynamical range 
of liquid relaxations. The deviations span about three 
orders of magnitude, compared with the overall range 
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FIG. 13: From Ref.[60]: Predictions for (T c - T K )/T K based 
on string-like rearranging regions is a dashed line, and (T c — 
Tk) /Tk based on percolation clusters is the solid line. The ex- 
perimentally derived crossover temperatures, (T^ xp —Tk) /Tk , 
from those materials collected by Novikov and Sokolov [63], 
are shown as circles with the dark circles referring to poly- 
mers. 



of 10 17 or so, accessible in experiment. At these faster 
times, according to Ref. [60], reconfigurations involving 
non-compact regions become important. This is consis- 
tent with "string-like" excitations that had been observed 
in numerical simulations of liquid dynamics by Glotzer 
and coworkers [64]. 

Generally simulations are carried out in the dynamic 
range where these softening effects are important. Yet, as 
in the last example, the RFOT theory is indeed consistent 
with current simulational studies. To be specific, cur- 
rent computational technologies are limited to a dynamic 
range of nine decades or so. The time step in simulation 
is roughly a femtosecond, significantly less than a typical 
collision time. In view of the large system size neces- 
sary to avoid boundary effects, at best one produces cur- 
rently a microsecond run, implying only thermally repre- 
sentative relaxations on a small fraction of a microsecond 
time-scale can be ascertained. (Equilibrating a liquid in 
a realistic simulation, in the deeply supercooled regime, 
must of course be subject to the very same high barriers 
present in the laboratory melts.) According to a recent 
discussion of Lubchenko [65], such short times imply the 
size of the critical nucleus (or the participation number at 
the saddle point, if you will) is only about (47r/3)l 3 ~ 5 
"beads", corresponding to r'/a ~ 1. In silica, the bead 
is a fraction of a SiC>4 tetrahedron [41, 66], implying a 
critical region size of 15 atoms or so. This is consistent 
with the participation ratios, at the barrier top, reported 
recently by Reinisch and Heuer in simulations character- 
istic of these short times [67]. According to this discus- 
sion, one must be careful so as not to mistake a possible 
but unrepresentative high barrier trajectory, that could 



be always found given a high enough temperature, for 
a thermally representative trajectory that actually de- 
scribes reconfiguration on the laboratory time scale. 



E. Aging 

A liquid is in equilibrium. Ultimately there is no mem- 
ory of its initially prepared configuration. But when the 
liquid is cooled faster than it can equilibrate, the sys- 
tem will find itself in a subset of all the states avail- 
able before the quench began. Which subset depends 
on the thermal history. The resultant quenched glass is 
a truly non-ergodic system. It will begin to relax to- 
wards the thermodynamic state which would have been 
typical upon slower quenching. But unlike at tempera- 
tures above vitrification, the structural re-arrangement 
below T g no longer depends only on the ambient tem- 
perature T alone but also on the temperature history. 
The vibrational temperature is near ambient, while the 
(now mostly static) structure is representative of an equi- 
librated sample at a "fictive" temperature Tf. Strictly, 
one must detail the complete temperature schedule. Yet, 
since the relaxations near a glass transition on the rou- 
tine laboratory scales are so slow - seconds and slower 
- that most quenches would result in a structure typi- 
cal of a single temperature T g , save the smaller vibra- 
tional amplitudes. This immediately implies the energy 
barrier distribution should be nearly temperature inde- 
pendent. If one further replaces the full barrier distri- 
bution by a single typical barrier, one expects the tem- 
perature dendence of relaxation times in frozen glasses 
to follow a simple Arrhenius law. A similar expectation 
serves as the basis of the phenomenological framework of 
Nayaranaswany-Moynihan-Tool: 

i, u J AE * n \ AE *\ 

K.e. = fc exp i 2? JVMT - — - (1 - X NM T)-r^T > ■ 

where E* is the equilibrated apparent activation energy 
at T g and inmt lies between and 1. The general ex- 
pression for activated transitions within RFOT theory 
in Eq.(14), which allows one to treat transitions from 
an arbitrary initial local state and implies the apparent 
activation energy below T g is significantly smaller than 
above the glass transition. The free energy difference 
fin ~ feq becomes nearly temperature independent be- 
low T g , also consistent with the initial state energy being 
above the equilibrated value. On the other hand, the rate 
of the temperature change of the apparent activation en- 
ergy above T g depends on the temperature dependence 
of the configurational entropy, see Eq.(16). The RFOT 
theory thus predicts the nonlinearity parameter xnmt to 
be correlated with the fragility. A straightforward calcu- 
lation yields the following simple relation [37] : 
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FIG. 14: The fragility parameter m is plotted as a function 
of the NMT nonlinearity parameter £nmt- The curve is pre- 
dicted by the RFOT theory when the temperature variation 
of 70 is neglected. The data are taken from Ref. [68]. The 
more fragile substances consistently lie above the prediction, 
which has no adjustable parameters. This discrepancy may 
be due to softening effects. 



This prediction is compared to experimental data in Fig- 
ure 14. (Here, we set Tj = T g .) In addition to predict- 
ing the typical activation energy in frozen glasses, the 
RFOT theory can be employed to estimate the spread 
in the barrier distribution, and thus enable one to com- 
pute the temperature dependence of the exponent (3 of 
the stretched exponential relaxation, as in Eq.(24). Cor- 
responding experimental data addressing this prediction 
so far appear inconclusive [69, 70]. 



IV. QUANTUM THEORY OF GLASSES 



implies there are local structural excitations that may be 
viewed as local multi-level systems. These anharmonic 
motions would interact with the mostly elastic lattice. At 
low enough temperatures, one expects only the two low- 
est energy levels to contribute thermodynamically (but 
see below). This is in line with the early, phcnomeno- 
logical two-level system (TLS) model which empirically 
describes the low T regime [73, 74]. But the low energy 
tail of the structural excitations can be computed from 
first principles. First one finds it is nearly energy inde- 
pendent at the lowest energies, as in the phenomenology. 
The corresponding density of states P is given by a sim- 
ple expression [38]: 



(29) 



This gives a linear heat capacity whose magnitude is con- 
sistent with experiment. The condition of marginal sta- 
bility of structural excitations to lattice distortions, at 
the glass transition, allows one to derive the coupling of 
those frozen-in transitions to the phonons: 



k B T g pc 2 a 3 , 



(30) 



where p and c 2 are the mass density and the speed of 
sound of the glass. These two expressions, represent- 
ing the nature of the states frozen at T g , enable us to 
understand the mysterious universality of phonon scat- 
tering, found in many amorphous insulators: The ratio 
of the phonon mean free path to the phonon wave length 
depends on Pg 2 and turns out to be about 150 for all 
substances. Thus the standard tunneling model implies 
Pg 2 /pc 2 s ~ 10~ 2 [75]. Combining Eqs.(29) and (30) from 
RFOT theory shows that the universality is expected be- 
cause it reflects the nearly universal value of the cooper- 
ativity length at the glass transition [38] : 



The microscopic picture of structural relaxations in su- 
percooled liquids and glasses built up from the RFOT 
theory provides a starting point for understanding a rad- 
ically different dynamical regime, i.e. the behavior of 
structural glasses at cryogenic temperatures [38], [39], 
[66], [71]. As we have recently provided a detailed 
account of these findings elsewhere [66, 71], we only 
summarize here several important points which help to 
strengthen the case for the RFOT theory predictions near 

T 
1 g- 

Much below the Debye temperature To, cryogenic 
glasses exhibit a set of excitations in excess of the usual 
Debye density of states. At temperatures below IK or so, 
these lead to a linear heat capacity and enhanced phonon 
scattering. In a higher temperature range around 10 K 
or so (~ 1 THz), yet another pronounced feature in the 
apparent density of excitations, called the "Boson Peak" , 
is observed [72] . This excess density of states also rapidly 
shortens the mean free path of thermal phonons. 

The mosaic predicted by RFOT theory immediately 



Pg 2 / P c 2 s 



1/200, 



(31) 



Another simple correlation empirically found but other- 
wise unexpected follows from RFOT theory. It relates 
the TLS-phonon coupling g to the glass transition tem- 
perature T g via the Lindcmann ratio: 



9 - 



(d L /a) 



~ 10.k B T„. 



(32) 



At the very lowest temperatures, the structural transi- 
tions involve the tunneling between two states of domain 
wall traversing a local region, spanning a length £ across. 
But the RFOT theory suggests these motions can sup- 
port more than two states. At temperatures above 1 K 
or so, higher energy states of these tunneling centers are 
predicted to come into play. Microscopically, these mo- 
tions can be visualized as vibrational excitations of the 
domain walls. The spectrum of these "capillary waves", 
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or "ripplons" , can be computed without adjustable pa- 
rameters. The corresponding frequencies and density of 
states quantitatively account for the apparent excess heat 
capacity and phonon scattering in the terahertz range, 
that have been associated with the Boson Peak. The 
simplest quantitative estimate of the vibration frequen- 
cies is obtained neglecting interaction with the phonons 
and yields for the Boson Peak frequency [39]: 

wbp — (a/£)<JD, (33) 

where wd is the Debye frequency. Effects of ripplon- 
phonon interaction will lead to a shift and brodening of 
the ripplon resonances, and seem to explain the non- 
universality of the thermal conductivity plateau. The 
multilevel nature of "two level systems" in this tem- 
perature range is also apparent in single molecule spec- 
troscopy experiments when complex spectral trails have 
ben observed by Orrit and coworkers [76]. 

The tunneling centers (TC) are fluctuating resonances 
and hence will mutually attract, following the usual Lon- 
don dispersion law r~ 6 . Since the number of active cen- 
ters grows with temperature, this implies the inter-center 
attraction provides a mechanism for an additional at- 
tracting holding the solid together at higher tempera- 
tures. Like rubber, this yields a negative thermal expan- 
sion coefficient! Quantitative estimates [66] show that the 
attraction is greatly enhanced by the presence of ripplons, 
even at very low T. This mechanism yields a magnitude 
of the Griineisen parameter consistent with experiment, 
in a number of glasses, which is often very anomalously 
large and negative [77]. 

Finally, individual molecular motions during the struc- 
tural transitions will result in local electric charge redis- 
tribution. Ultimately this leads to a coupling of the tran- 
sitions with electromagnetic waves. Recent estimates [71] 
show the transition induced electric dipole is about a De- 
bye in magnitude, in spite of there being several hundred 
atoms participating in the motion. This value is consis- 
tent with experiments [78, 79]. 

V. CONCLUDING REMARKS - FUTURE 
DIRECTIONS 

In place of a summary, we discuss the range of appli- 
cability of the RFOT theory; the limitations uncovered 
naturally suggest directions for future efforts. 

We have already pointed out that the RFOT theory's 
predictions of the relaxation time and size scales are 
quantitative at viscosities of 10 Poise and above. This 
range of time scales constitutes, on the logarithmic scale, 
at least 80% of the routinely probed dynamic range in 
liquids of 10~ 2 to 10 14 Poise or so. The ability of the 
theory to make predictions here relies on the fact that 
the motions are mesoscopic, hundreds of atoms being re- 
configured, allowing sufficient thermodynamic averaging. 

Based on mesoscale averaging the theory thus applies 
only when there is no significant medium range order. 



Many liquids seem to satisfy this constraint, and span 
a wide range fragilities and chemical properties. These 
include nearly ionic compounds, such SiC>2, ZnCl2; al- 
cohols, such as propanol; aromatic and aliphatic hydro- 
carbons such as trinaphthylbenzene or 2-mcthylpcntanc. 
Measurements on such simple liquids confirm the predic- 
tions of the RFOT theory. The list of substances quan- 
titatively described presently includes several dozens of 
materials and keeps growing. 

On the other hand, some substances, while showing 
most of the signatures of glassy behaviors, exhibit addi- 
tional types of order and scales of motion that will require 
separate treatment. We have in fact mentioned several 
examples of such systems: They include many of the most 
important polymers. In polymers, for instance, it is clear 
that various chain motions may provide long length scales 
and long time-scale relaxations, especially when the per- 
sistence length is large. Even for such systems, one still 
finds qualitative agreement with the RFOT calculations. 
For example, most polymers show TLS and Boson Peak 
densities of states that are comparable to those of simple 
molecular systems. One consequence of the possiblity 
of mid-range order is that one should exercise caution 
in searching for universalities in the glassy behaviours. 
The RFOT theory, among other things, provides in fact 
a first principles basis for such searches. As we have seen, 
substances that conform to the RFOT specifications, are 
predicted to show not one but rather a whole set of corre- 
lations, most of which we have mentioned earlier: Tk vs. 
To, Ac p vs. D (or to), £ vs. T, £(T g ) vs. a, universality 
of s c (T g )_, (3 vs. D, 13 vs. T, m vs. x, T c vs. Ac p , P 
vs. Tg, P vs. £ (or a), P vs. g, cjbp vs. ivd, ^bp vs. a. 
The experiments become more difficult as one proceeds 
down the list, so fewer comparisons of predictions with 
experiment are available for the last few relations, but in 
our view the number of confirmed relationships already 
provides a secure basis to suggest anomalies are a sign of 
new degrees of freedom. 

We were careful to limit the title of this review calling it 
"Theory of Structural Glasses and Supercooled Liquids" , 
not dealing with all glassy systems. The quantitative 
applicability of RFOT theory to molecular glassy sys- 
tems can be ultimately traced back to an emergent small 
parameter, the Lindcmann ratio. Also for structural 
glasses, the glass transition temperature T g is much be- 
low the dynamic transition at Ta- This was pointed out 
clearly in an analysis by Eastwood and Wolyncs [80] who 
developed a "Ginzburg" criterion for the theory. Their 
criterion shows many exciting systems exhibiting signs 
of glassiness are not in the appropriate universal regimes 
for strict quantitation because those systems are intrinsi- 
cally softer with larger effective local motions. Colloidal 
glasses, owing to the larger size of their constituents, are 
intrinsically slow at the one particle level and thus at hu- 
man measurement time scales, are always near Ta- Gels, 
stripe glasses [30], and microemulsions [81] have larger 
Lindcmann ratios and although they are described by 
the RFOT theory at the mean-field level, again, we ex- 
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pect there will be very significant renormalization of their 
behavior in the activated regime. Finally, although the 
RFOT theory provides a route from input intermolccular 
forces to the macroscopic behavior, the present quantita- 
tive successes of the theory avoid the hardest part of the 
detailed microscopic modeling by utilizing the fact that it 
is relatively easy to estimate conhgurational entropy from 
laboratory measurements. As usual, completely ab initio 



calculations will be more difficult, justifying greater at- 
tention to microscopic liquid state theory and molecular 
simulation technology. 
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